Diagnosing destabilization risk in global land carbon sinks.

Global net land carbon uptake or net biome production (NBP) has increased during recent decades1. Whether its temporal variability and autocorrelation have changed during this period, however, remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink2,3. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO2 concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.

Future bioenergy expansion could alter carbon sequestration potential and exacerbate water stress in the United States.

The maximum future projected bioenergy expansion potential, in scenarios limiting warming to 2°C or below, is equivalent to half of present-day croplands. We quantify the impacts of large-scale bioenergy expansion against re/afforestation, which remain elusive, using an integrated human-natural system modeling framework with explicit representation of perennial bioenergy crops. The end-of-century net carbon sequestration due to bioenergy deployment coupled with carbon capture and storage largely depends on fossil fuel displacement types, ranging from 11.4 to 31.2 PgC over the conterminous United States. These net carbon sequestration benefits are inclusive of a 10 PgC carbon release due to land use conversions and a 2.4 PgC loss of additional carbon sink capacity associated with bioenergy-driven deforestation. Moreover, nearly one-fourth of U.S. land areas will suffer severe water stress by 2100 due to either reduced availability or deteriorated quality. These broader impacts of bioenergy expansion should be weighed against the costs and benefits of re/afforestation-based strategies.

Increasing the spatial and temporal impact of ecological research: A roadmap for integrating a novel terrestrial process into an Earth system model.

Terrestrial ecosystems regulate Earth's climate through water, energy, and biogeochemical transformations. Despite a key role in regulating the Earth system, terrestrial ecology has historically been underrepresented in the Earth system models (ESMs) that are used to understand and project global environmental change. Ecology and Earth system modeling must be integrated for scientists to fully comprehend the role of ecological systems in driving and responding to global change. Ecological insights can improve ESM realism and reduce process uncertainty, while ESMs offer ecologists an opportunity to broadly test ecological theory and increase the impact of their work by scaling concepts through time and space. Despite this mutualism, meaningfully integrating the two remains a persistent challenge, in part because of logistical obstacles in translating processes into mathematical formulas and identifying ways to integrate new theories and code into large, complex model structures. To help overcome this interdisciplinary challenge, we present a framework consisting of a series of interconnected stages for integrating a new ecological process or insight into an ESM. First, we highlight the multiple ways that ecological observations and modeling iteratively strengthen one another, dispelling the illusion that the ecologist's role ends with initial provision of data. Second, we show that many valuable insights, products, and theoretical developments are produced through sustained interdisciplinary collaborations between empiricists and modelers, regardless of eventual inclusion of a process in an ESM. Finally, we provide concrete actions and resources to facilitate learning and collaboration at every stage of data-model integration. This framework will create synergies that will transform our understanding of ecology within the Earth system, ultimately improving our understanding of global environmental change, and broadening the impact of ecological research.

Interannual and Seasonal Drivers of Carbon Cycle Variability Represented by the Community Earth System Model (CESM2).

Earth system models are intended to make long-term projections, but they can be evaluated at interannual and seasonal time scales. Although the Community Earth System Model (CESM2) showed improvements in a number of terrestrial carbon cycle benchmarks, relative to its predecessor, our analysis suggests that the interannual variability (IAV) in net terrestrial carbon fluxes did not show similar improvements. The model simulated low IAV of net ecosystem production (NEP), resulting in a weaker than observed sensitivity of the carbon cycle to climate variability. Low IAV in net fluxes likely resulted from low variability in gross primary productivity (GPP)-especially in the tropics-and a high covariation between GPP and ecosystem respiration. Although lower than observed, the IAV of NEP had significant climate sensitivities, with positive NEP anomalies associated with warmer and drier conditions in high latitudes, and with wetter and cooler conditions in mid and low latitudes. We identified two dominant modes of seasonal variability in carbon cycle flux anomalies in our fully coupled CESM2 simulations that are characterized by seasonal amplification and redistribution of ecosystem fluxes. Seasonal amplification of net and gross carbon fluxes showed climate sensitivities mirroring those of annual fluxes. Seasonal redistribution of carbon fluxes is initiated by springtime temperature anomalies, but subsequently negative feedbacks in soil moisture during the summer and fall result in net annual carbon losses from land. These modes of variability are also seen in satellite proxies of GPP, suggesting that CESM2 appropriately represents regional sensitivities of photosynthesis to climate variability on seasonal time scales.

Towards a multiscale crop modelling framework for climate change adaptation assessment.

Predicting the consequences of manipulating genotype (G) and agronomic management (M) on agricultural ecosystem performances under future environmental (E) conditions remains a challenge. Crop modelling has the potential to enable society to assess the efficacy of G × M technologies to mitigate and adapt crop production systems to climate change. Despite recent achievements, dedicated research to develop and improve modelling capabilities from gene to global scales is needed to provide guidance on designing G × M adaptation strategies with full consideration of their impacts on both crop productivity and ecosystem sustainability under varying climatic conditions. Opportunities to advance the multiscale crop modelling framework include representing crop genetic traits, interfacing crop models with large-scale models, improving the representation of physiological responses to climate change and management practices, closing data gaps and harnessing multisource data to improve model predictability and enable identification of emergent relationships. A fundamental challenge in multiscale prediction is the balance between process details required to assess the intervention and predictability of the system at the scales feasible to measure the impact. An advanced multiscale crop modelling framework will enable a gene-to-farm design of resilient and sustainable crop production systems under a changing climate at regional-to-global scales.

Challenging a Global Land Surface Model in a Local Socio-Environmental System.

Land surface models (LSMs) predict how terrestrial fluxes of carbon, water, and energy change with abiotic drivers to inform the other components of Earth system models. Here, we focus on a single human-dominated watershed in southwestern Michigan, USA. We compare multiple processes in a commonly used LSM, the Community Land Model (CLM), to observational data at the single grid cell scale. For model inputs, we show correlations (Pearson's R) ranging from 0.46 to 0.81 for annual temperature and precipitation, but a substantial mismatch between land cover distributions and their changes over time, with CLM correctly representing total agricultural area, but assuming large areas of natural grasslands where forests grow in reality. For CLM processes (outputs), seasonal changes in leaf area index (LAI; phenology) do not track satellite estimates well, and peak LAI in CLM is nearly double the satellite record (5.1 versus 2.8). Estimates of greenness and productivity, however, are more similar between CLM and observations. Summer soil moisture tracks in timing but not magnitude. Land surface reflectance (albedo) shows significant positive correlations in the winter, but not in the summer. Looking forward, key areas for model improvement include land cover distribution estimates, phenology algorithms, summertime radiative transfer modelling, and plant stress responses.

Dry Deposition of Ozone over Land: Processes, Measurement, and Modeling.

Dry deposition of ozone is an important sink of ozone in near surface air. When dry deposition occurs through plant stomata, ozone can injure the plant, altering water and carbon cycling and reducing crop yields. Quantifying both stomatal and nonstomatal uptake accurately is relevant for understanding ozone's impact on human health as an air pollutant and on climate as a potent short-lived greenhouse gas and primary control on the removal of several reactive greenhouse gases and air pollutants. Robust ozone dry deposition estimates require knowledge of the relative importance of individual deposition pathways, but spatiotemporal variability in nonstomatal deposition is poorly understood. Here we integrate understanding of ozone deposition processes by synthesizing research from fields such as atmospheric chemistry, ecology, and meteorology. We critically review methods for measurements and modeling, highlighting the empiricism that underpins modeling and thus the interpretation of observations. Our unprecedented synthesis of knowledge on deposition pathways, particularly soil and leaf cuticles, reveals process understanding not yet included in widely-used models. If coordinated with short-term field intensives, laboratory studies, and mechanistic modeling, measurements from a few long-term sites would bridge the molecular to ecosystem scales necessary to establish the relative importance of individual deposition pathways and the extent to which they vary in space and time. Our recommended approaches seek to close knowledge gaps that currently limit quantifying the impact of ozone dry deposition on air quality, ecosystems, and climate.

Beyond Static Benchmarking: Using Experimental Manipulations to Evaluate Land Model Assumptions.

Land models are often used to simulate terrestrial responses to future environmental changes, but these models are not commonly evaluated with data from experimental manipulations. Results from experimental manipulations can identify and evaluate model assumptions that are consistent with appropriate ecosystem responses to future environmental change. We conducted simulations using three coupled carbon-nitrogen versions of the Community Land Model (CLM, versions 4, 4.5, and-the newly developed-5), and compared the simulated response to nitrogen (N) and atmospheric carbon dioxide (CO2) enrichment with meta-analyses of observations from similar experimental manipulations. In control simulations, successive versions of CLM showed a poleward increase in gross primary productivity and an overall bias reduction, compared to FLUXNET-MTE observations. Simulations with N and CO2 enrichment demonstrate that CLM transitioned from a model that exhibited strong nitrogen limitation of the terrestrial carbon cycle (CLM4) to a model that showed greater responsiveness to elevated concentrations of CO2 in the atmosphere (CLM5). Overall, CLM5 simulations showed better agreement with observed ecosystem responses to experimental N and CO2 enrichment than previous versions of the model. These simulations also exposed shortcomings in structural assumptions and parameterizations. Specifically, no version of CLM captures changes in plant physiology, allocation, and nutrient uptake that are likely important aspects of terrestrial ecosystems' responses to environmental change. These highlight priority areas that should be addressed in future model developments. Moving forward, incorporating results from experimental manipulations into model benchmarking tools that are used to evaluate model performance will help increase confidence in terrestrial carbon cycle projections.

Comparing optimal and empirical stomatal conductance models for application in Earth system models.

Earth system models (ESMs) rely on the calculation of canopy conductance in land surface models (LSMs) to quantify the partitioning of land surface energy, water, and CO2 fluxes. This is achieved by scaling stomatal conductance, gw , determined from physiological models developed for leaves. Traditionally, models for gw have been semi-empirical, combining physiological functions with empirically determined calibration constants. More recently, optimization theory has been applied to model gw in LSMs under the premise that it has a stronger grounding in physiological theory and might ultimately lead to improved predictive accuracy. However, this premise has not been thoroughly tested. Using original field data from contrasting forest systems, we compare a widely used empirical type and a more recently developed optimization-type gw model, termed BB and MED, respectively. Overall, we find no difference between the two models when used to simulate gw from photosynthesis data, or leaf gas exchange from a coupled photosynthesis-conductance model, or gross primary productivity and evapotranspiration for a FLUXNET tower site with the CLM5 community LSM. Field measurements reveal that the key fitted parameters for BB and MED, g1B and g1M, exhibit strong species specificity in magnitude and sensitivity to CO2 , and CLM5 simulations reveal that failure to include this sensitivity can result in significant overestimates of evapotranspiration for high-CO2 scenarios. Further, we show that g1B and g1M can be determined from mean ci /ca (ratio of leaf intercellular to ambient CO2 concentration). Applying this relationship with ci /ca values derived from a leaf δ13 C database, we obtain a global distribution of g1B and g1M , and these values correlate significantly with mean annual precipitation. This provides a new methodology for global parameterization of the BB and MED models in LSMs, tied directly to leaf physiology but unconstrained by spatial boundaries separating designated biomes or plant functional types.

Monoterpene 'thermometer' of tropical forest-atmosphere response to climate warming.

Tropical forests absorb large amounts of atmospheric CO2 through photosynthesis but elevated temperatures suppress this absorption and promote monoterpene emissions. Using 13 CO2 labeling, here we show that monoterpene emissions from tropical leaves derive from recent photosynthesis and demonstrate distinct temperature optima for five groups (Groups 1-5), potentially corresponding to different enzymatic temperature-dependent reaction mechanisms within ß-ocimene synthases. As diurnal and seasonal leaf temperatures increased during the Amazonian 2015 El Niño event, leaf and landscape monoterpene emissions showed strong linear enrichments of ß-ocimenes (+4.4% °C-1 ) at the expense of other monoterpene isomers. The observed inverse temperature response of α-pinene (-0.8% °C-1 ), typically assumed to be the dominant monoterpene with moderate reactivity, was not accurately simulated by current global emission models. Given that ß-ocimenes are highly reactive with respect to both atmospheric and biological oxidants, the results suggest that highly reactive ß-ocimenes may play important roles in the thermotolerance of photosynthesis by functioning as effective antioxidants within plants and as efficient atmospheric precursors of secondary organic aerosols. Thus, monoterpene composition may represent a new sensitive 'thermometer' of leaf oxidative stress and atmospheric reactivity, and therefore a new tool in future studies of warming impacts on tropical biosphere-atmosphere carbon-cycle feedbacks.